Insulating spin liquid in the lightly doped two-dimensional Hubbard model 
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We calculate the charge compressibility and uniform spin susceptibility for the two-dimensional 
(2D) Hubbard model slightly away from half-filling within a two-loop renormalization group scheme. 
We find numerically that both those quantities flow to zero as we increase the initial interaction 
strength from weak to intermediate couplings. This result implies gap openings in both charge and 
spin excitation spectra for the latter interaction regime. When this occurs, the ground state of the 
lightly doped 2D Hubbard model may be interpreted as an insulating spin liquid as opposed to a 
Mott insulating state. 
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After two decades of intensive research on the high- 
Tc superconductors, physicists are still puzzled by some 
of their very unusual electronic properties^. The promi- 
nent example is given by the cuprates. At zero doping, 
despite the fact that their highest occupied band is half- 
filled, they are charge insulators, and display antiferro- 
magnetic long-range order. For this reason, they are said 
to be Mott insulators. As soon as one starts doping those 
compounds with holes, the long-range magnetic order be- 
comes rapidly suppressed, and there are experimental ev- 
idences of an emergence of a spin gap in their correspond- 
ing excitation spectra™. A charge gap is also observed 
by ARPES experiments in such lightly-doped systems*. 
Moreover, at finite temperatures, they turn themselves 
into poor conductors with electronic properties differing 
considerably from the predictions of Landau's Fermi liq- 
uid theory. This scenario configures the so-called pseu- 
dogap regime. Although this phase continues to be not 
well understood, it is widely acknowledged to play a fun- 
damental role in the underlying microscopic mechanism 
of such high-Tc superconductors. Indeed, upon some fur- 
ther doping, those poor metals become superconducting 
with d-wave symmetry up to relatively high temperatures 
around the optimal doping level. 

From the theoretical viewpoint, it is widely accepted 
that the appropriate model for describing such systems is 
the two-dimensional (2D) Hubbard model (HM), since it 
is known to have a Mott insulating phase at half-filling, 
and is expected to become a d-wave superconductor at 
larger doping 4 . However, its intermediate doping regime, 
which could provide some insight to understand the phys- 
ical nature of the pseudogap state, still remains elusive 
to this date. In this Letter we intend to address this 
question using renormalization group (RG) techniques in 
order to infer about the ground state of such model for 
electron densities slightly away from the half- filling limit. 

Our considerations here will be based on a complete 
two-loop RG calculation of the uniform charge (CS) and 
spin (SS) susceptibilities of the 2D HM, taking into ac- 
count simultaneously both the renormalization of the 
couplings, and the self-energy effects. (The CS is also 
called the charge compressibility of the system.) To 
best of our knowledge, it is the first time that such a 



full two-loop RG calculation is performed for the 2D 
lightly-doped HM, since previous estimates of the uni- 
form susceptibilities followed random phase approxima- 
tion (RPA) schemes. 

In momentum space, the 2D Hubbard Hamiltonian on 
a square lattice is given by 

k.cr ^ si es / p j^ q 
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where the energy dispersion is simply £k = 
— 2t [cos^a) + cos(k y a)] — /i, and ip^ a and i/w 
are the usual creation and annihilation operators of 
electrons with momentum k and spin projection er =f, J,. 
Besides, /i stands for the chemical potential, whereas 
a is the square lattice spacing. Another important 
parameter here is the width of the noninteracting band, 
which is given by W = 8t. This model describes a 
system with many electrons interacting mutually via a 
local repulsive interaction U, and with a total number 
Nsites of lattice sites. The electron band filling of the 
system is controlled by the ratio fi/t. When (x/t = the 
system is exactly at half-filling. As we start doping it 
with holes, /i/i takes slightly negative values. 
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FIG. 1: (Color online) (a) The half-filled Fermi surface (FS) of 
the 2D Hubbard model (dashed line), the lightly hole-doped 
FS (blue solid line), and (b) the latter FS after a rotation of 
the axes by 45° degrees. 



Our starting point is a 2D nearly flat Fermi surface 
(FS) with no van Hove singularities (see Fig. 1). It cor- 
rectly describes the 2D HM slightly away from the half- 
filling case. A similar FS has already been used by other 
groups to investigate the leading instabilities within ei- 
ther a parquet or, equivalently, a one-loop RG approach 5 . 
All those investigations find diverging susceptibilities at 
finite energies (or finite temperatures) with the dominant 
instability being always the spin density wave (SDW). 
Their interpretation is that this implies a spontaneous 
symmetry breaking in the system, and the onset of a 
long-range ordered antiferromagnetic state. 

In this Letter we argue that is not necessarily true 
since this result may also be an indicative of the limi- 
tations of the one-loop RG scheme. In low-dimensional 
systems, large quantum fluctuations are expected to sup- 
press long-range order. The more those effects are taken 
into account, the more likely those long-range ordered 
states are transformed into short-range magnetically or- 
dered phases. This also becomes clear as a result of 
our work. Taking into account quantum fluctuation ef- 
fects up to two-loop order in our RG scheme, we are 
able to show that, for moderate coupling regimes, both 
charge compressibility and uniform spin susceptibility are 
strongly suppressed and flow unequivocally to zero. This 
behavior implies that there are gaps for both charge and 
spin excitations, and no trace of long-range symmetry 
breaking order in those cases. Such a state with a fully 
gapped charge and spin spectra is usually denominated 
an insulating spin liquid (ISL). The ISL is an example 
of a short-range resonant valence bond state, which was 
first proposed for a S=l/2 Heisenberg mode!-., and clearly 
revealed in even-leg Hubbard ladders by both RG and 
bosonization approaches'^. 

In order to implement a full RG calculation of the uni- 
form susceptibilities, it is essential to consider at least 
two-loop order contributions. This is due to the fact that, 
at one-loop level, there is not a single infrared (IR) diver- 
gent diagram in the calculation of the so-called uniform 
response functions. In contrast, there are several of those 
diagrams in two loops (the nonparquet diagrams), and, 
as a result, one can reliably begin to derive appropriate 
RG flow equations for those quantities at that order of 
perturbation theory. 

To keep a closer contact with well-known works in one- 
dimensional systems^ ", we divide the FS into four differ- 
ent regions (two sets of solid and dashed line patches). 
Here we restrict the momenta at the FS to the flat 
parts only. The interaction processes connecting paral- 
lel patches of the FS are always logarithmically IR di- 
vergent due to quantum fluctuations. In contrast, those 
connecting perpendicular patches always remain finite, 
and do not contribute to the RG flow equations in our 
approach. For convenience, we restrict ourselves to one- 
electron states labeled by the momenta pn = k x and 
p± = k y associated with one of the two sets of per- 
pendicular patches. The momenta parallel to the FS 
are restricted to the interval —A ^ p» sC A, with 2A 
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FIG. 2: The uniform response function r M (p,q«0) up to 
two-loop order. The triangular vertices stand for the renor- 
malized response function, whereas the diagram with a cross 
represents the counterterm. 

being essentially the size of the flat patches. The en- 
ergy dispersion of the single-particle states is given by 
£ a (p) = vp ([px| — hp), and depends only on the mo- 
menta perpendicular to the FS. The label a = ± refers 
to the flat sectors at p± = ±fc^ , respectively. In addi- 
tion, we take kp — A ^ \p±\ ^ kp + A, where A is a fixed 
ultraviolet (UV) microscopic momentum cut-off. 

We now write down the Lagrangian of the 2D HM as 

L= H (P) ~Vf(\p±\ -k F )]ip (a)(r (p) 

p,cr,a— ± 

p,q,k a,0,5,-y 

x ^l +) s (p + q - k ) ^(- )7 O) i>(.-)p (q) ^(+)a (p) » 

(2) 

where W±) an d ^(±) are now fermionic fields associ- 
ated to electrons located at the ± patches. The sum- 
mation over momenta must be appropriately understood 
as J2 P — V/i^) 2 J d 2 p in the thermodynamic limit. 
We linearized the energy dispersion about the lightly- 
doped FS, and the interaction term was parametrized 
in a manifestly SU(2) invariant form. Here we follow 
the well-known g-ology notation, with g\ and <?2 stand- 
ing for backscattering and forward scattering couplings, 
respectively. Since this should represent the 2D HM, 
these couplings must be initially defined as 31 = 32 = 
(V '/ 'AN S ites)U . In addition, we do not include Umklapp 



3 



processes since we are not at half-filling, and our FS does 
not intersect the so-called Umklapp surface at any point 
(see Fig. 1(a)). Consequently, our result is different from 
another evidence of ISL behavior reported in the litera- 
ture in the context of the t — t' 2D HM 11 . In their case, 
the FS intersects the Umklapp surface even away from 
half-filling, and, for this reason, Umklapp processes be- 
come an essential ingredient for the correct description 
of the system at the lightly-doped regime. Moreover, 
their estimates of the uniform susceptibilities follow the 
already mentioned RPA-like scheme. 

Since the HM is a microscopic model, all terms in the 
Lagrangian are defined at a scale of a few lattice spac- 
ings in real space (i.e at the UV cutoff scale in momen- 
tum space). These parameters are often inaccessible to 
every day experiments, since the latter probe only the 
low-energy and the long-wavelength dynamics of a given 
system. In RG theory, these unobserved quantities are 
known as the bare parameters. In fact, if one attempts to 
construct naive perturbative calculations with such pa- 
rameters, one will obtain IR divergent Feynman diagrams 
in the computation of several quantities such as, e.g. the 
backscattering and the forward scattering four-point ver- 
tices, and the single-particle Green's function. These di- 
vergences mean that the perturbation theory setup is not 
appropriately formulated. To solve this problem, one can 
redefine the perturbation scheme in a such a way as to 
circumvent these infinite results in the calculation of ob- 
servable quantities. This is the strategy of the so-called 
renormalized perturbation theory. Thus, one rewrites all 
unobserved bare parameters in terms of the associated 
observable renormalized (or physical) quantities. The 
difference between them is given by a counterterm, whose 
fundamental role is to cancel, by construction, the related 
IR divergences in all orders of perturbation theory. If this 
program is successfully accomplished, then the theory is 
said to be properly renormalized. 

As we explained in our previous paper—, when one 
deals with this FS problem, the counterterms needed to 
renormalize the theory turn out to be continuous func- 
tions of the three momenta parallel to the FS, rather 
than being simply infinite constants. In addition, we 
computed the RG flow equations for the coupling func- 
tions and the quasiparticle weight up to two-loop or- 
der, and showed that they were in fact coupled integro- 
differential equations. We solved those equations self- 
consistently and, as a result, we found out for an inter- 
mediate coupling regime a possibly new physical regime, 
which was characterized by a strongly suppressed quasi- 
particle weight. Here we continue to explore this inter- 



acting regime, and our present calculation shows that the 
resulting quantum state may in fact be interpreted as an 
ISL. 

To obtain the uniform susceptibilities of this system, 
we must first calculate the linear response due to an 
infinitesimal external field perturbation that couples to 
both charge and spin number operators. We do this by 
adding to the Lagrangian the new term 



-K xtmal Y, ^P)<] Q (p)Vf a)Q (p), (3) 

p,a— ± 

where B stands for the bare quantities. This will generate 
an additional vertex (the one-particle irreducible function 
r^ 2,1 -* (p, q w 0)), which will in turn be afflicted by new 
IR divergences (see the nonparquet diagrams in Fig. 2). 
As a result, we must rewrite the bare quantity Tg a in 
terms of its renormalized counterpart (henceforth called 
Tft a ), and an appropriate counterterm AT^ Q , i.e. 



Tg° On) = z-\ v {) [rr M ) + Arf (pi,)] . (4) 

The quasiparticle weight Z factor comes from the renor- 
malization of the fermionic fields, which must be also 
taken into account in order to include the feedback of 
the self-energy effects into the RG flow equations. The Z 
function is calculated explicitly in Ref— . As was men- 
tioned before, AT^ Q must cancel exactly the divergences 
generated by the nonparquet diagrams. However, there 
are still several ways of choosing that counterterm. To 
solve this ambiguity, we must make a prescription es- 
tablishing precisely that the T^ a is the experimentally 
observable response, i.e., r( 2,1 )(p||,po = u,p± = &F;q ~ 
0) = — iTft a (p\\,Lu), where u> is the RG energy scale pa- 
rameter that denotes the proximity of the renormalized 
theory to the FS. In this way, to flow towards the FS we 
let w -» 0. 

We are now ready to define the two different types 
of uniform response functions, which arise from a sym- 
metrization with respect to the spin projection, namely 

T R}C s(p h u) = T£(p h u) + TU(p h u), (5) 
T^ssip^uj) = 72 T (P||,w)-^(P||,w), (6) 
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FIG. 3: Feynman diagram associated with the renormalized 
charge compressibility and uniform spin susceptibility. 



where Tr.cs and Tr^ss are the response functions as- 
sociated with the charge compressibility and spin sus- 
ceptibility, respectively. In order to compute the RG 
flow equations for these response functions, one needs 
to recall that the bare parameters are independent of the 
renormalization scale u>. Thus, using the RG condition 
uidTg a jdto — 0, we obtain 
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FIG. 4: (Color online) RG flow of both charge compressibil- 
ity and uniform spin susceptibility as we increase the initial 
coupling strength given by U = (g/nvpjt. 

w4-Tk,cs(p\\) = -u-^AT r ,cs(p\\)+j{p\\)Tr,cs(p\\), 

(7) 

-Tj^ssipw) = -w— A7h,ss(p||) +j{p\\)T Ri ss{p\\), 

(8) 

where the anomalous dimension is given by ~f(p») = 
wd\nZ{p\\)/duJ. Despite their apparent simplicity, it is 
impossible to solve these RG equations only by analyti- 
cal means. To find their solutions, we have again to resort 
to numerics. Here we use the fourth-order Runge-Kutta 
numerical method. We discretize the FS continuum re- 
placing the interval —A ^ pu ^ A by a discrete set of 
33 points. We use to — £!exp(— /), where fl = 2vf^ with 
I being our RG step. We also choose f2 = vpA < W. 
In view of our choice of points for the FS, we are only 
allowed to go up to I w 2.8 in the RG flow to avoid the 
distance to the FS being smaller than the shortest dis- 
tance between points in our discrete set. 

Once the response functions are obtained, we can cal- 
culate the flow of the charge compressibility and the uni- 
form spin susceptibility of the system. Using again our 
diagrammatic convention, it follows from Fig. 3 that they 
are given by 



To evaluate these quantities we follow the same numeri- 
cal procedure described above. Our results are displayed 
in Fig. 4. We note in this plot that both charge and spin 
susceptibilities flow at the same rate as we approach the 
FS even though their corresponding response functions 
have different flow equations. In addition, for initial cou- 
plings in which the quasiparticle weight flows to zero, 
the uniform susceptibilities become strongly suppressed 
in the low-energy limit. While the former asserts that 
there are no fermionic quasiparticle excitations present 
in the system, the latter is related to the complete ab- 
sence of low- lying bosonic charge and spin collective exci- 
tations. Since this resulting state has only gapful excita- 
tions, it cannot be related to any broken symmetry state 
and, as a consequence, should possess only short-range 
ordering. This quantum state has, therefore, strong simi- 
larities with that predicted long ago by Anderson^, which 
is commonly referred to as an ISL. In our present work, 
it becomes evident that such a state is produced by dis- 
ordering effects induced by strong quantum fluctuations, 
and these are approximately taken into account by our 
two-loop RG scheme. Finally, we call attention to the 
fact that an insulating behavior in the lightly doped 2D 
HM was also reported recently in the literature^. Our 
present result is clearly in agreement with their results. 



In summary, we examined the flow of both charge com- 
pressibility and uniform spin susceptibility in the lightly 
doped 2D HM as a function of the initial interaction 
strength within a two-loop RG approach. For moderate 
interaction regimes, both quantities flow to zero as we 
approach the initial FS of the system. This is a strong 
indicative that there are gaps in both charge and spin 
excitation spectra of the lightly doped 2D HM. Hence 
the quantum state associated with that regime may be 
viewed as an ISL as opposed to a Mott insulator. This 
result may be of relevance for the cuprate high-Tc su- 
perconductors in view of the fact that the 2D HM in 
the intermediate coupling regime is widely believed to 
be appropriate to describe such compounds in all doping 
ranges. 
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